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ABSTRACT 

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type 
galaxies are close to isothermal, i.e., ptot °c r 1 with 7 w —2. To understand the origin of this universal 
slope we study a set of simulated spheroids formed in isolated binary mergers with controlled initial 
conditions as well as the formation within the cosmological framework. On average, the total stellar 
plus dark matter density profiles can be described by a power law with an index of 7 w —2.1 with 
a tendency towards steeper slopes for more compact, lower-mass ellipticals. In the binary mergers 
the amount of gas involved in the merger determines the steepness of the slope. This agrees with 
results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction 
of stars formed in-situ. At higher redshifts, the slopes of the ellipticals extracted from the cosmological 
simulations are generally steeper. Each gas-poor merger event evolves the slope towards 7 ~ —2, once 
this slope is reached further merger events do not change the slope anymore. Independent of their 
individual slopes or evolution scenarios, all our ellipticals have flat intrinsic combined stellar and dark 
matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density 
distributions with a slope of 7 ~ —2 for the combined system of stars and dark matter act as a 
natural attractor. In addition, the variety of complex formation histories as present in cosmological 
simulations, including major as well as minor merger events, is essential to generate the full range of 
observed density slopes seen for present day elliptical galaxies. 

Subject headings: galaxies: interactions - evolution - formation - cosmology: dark matter - methods: 
numerical 



1. INTRODUCTION 

Early-type galaxies are among the most massive 
and prominent galaxies of the universe, and thus 
their origin and formation history has been the fo- 
cus of many o bservatio n al st udies as well as simu- 
lations. Since iToomrd (|1977j ) proposed that ellipti- 
cal galaxies can be produced by a merger between 
two spiral galaxies, several simulations have st udied 
this mechanism in detail , for example IWhitd (I1978L 
1979aTEh: iGerhardl (fl98lh : iNeeroponte fc White! (fl 983): 



Barnes! (Il988f). including var i ous ga s fractions (|Hernquis 
19891: lBarnes_fc_Hernquistl 119961 : iNaab et all 12006 



Novak et al.ll2012l) and black h ole physics (Springcl et al.l 



2005al Uohansson et alj|2009al lbf) 



Although major merger simulations of spiral galax- 
ies could explain many observed properties of ellipti- 
cals, several problems were also identified. A number 
of studies have indicated that observations of massive, 
slowly rotating and spheroidal ellipticals cannot be ex- 
plained by wet or dry major mergers o f disk galaxies 
(INaab fc BurkertJ 120031: ICox et all 120061: iBurkert et all 
120081 : iBois et al.l 120 1 Ct 120 1 lh . Multiple early-type merg- 
ers would be needed in order to form the most mas- 
sive elliptica ls with masses that exceed the mass of sp i- 
ral galaxies llNaab et alJ!2006bL iKhochfar fc Silld f2006h . 
iGenel et al.l ( 201 Oh however demonstrated that the frac- 
tion of observed major early- type mergers is not high 
enough to explain the number of observed giant ellipti- 



cal galaxies. An additional complication is that the high 
metallicity and old ages observed in present day mas- 
sive early-type galaxies cannot be explained by binary 
mergers between tw o typical present day sp iral galaxies 
or their progenitors (|Naab fc Ost rikcr 200jJ see however 
iHopkins et aLll2008f ). 

Simulations of the formation of elliptical galaxies in 



all picture (e.g.. iMeza et al.l l2003t INaab et al 





2007 




2012 



Johansson et alJl2012t lLackner et al.ll2^12l ) , verifying the 
idea that elliptical galaxies can be produced by multiple 
minor mergers, a much more likel y evolution of events i n 
the lifetime of a galaxy (see also iBournaud et al.ll2007h . 
Mass accretion histories from cosmological simul ations 
(e.g.. IGao et"alll200l IFakhouri fc Mall2008l l200l have 
shown that the most likely process for forming elliptical 
galaxies is a mixture of both scenarios: most massive 
halos go through an early dissipative phase of fast accre- 
tion, including major merger events, followed by a phase 
of stellar accretion, during which all kind of mergers 
and even smooth accretion can be the dominant growth 
mechanism (e.g. JOser et al.ll2010D . 

There have been several approaches to connect the 
observed elliptical properties to the different forma- 
tion scenarios to investigate if we can explain observed 
features by the ev olution hist o ry of i ndividual ellipti- 
cals. For example. INaab et al.1 (jl999); B endo fc Barnes! 
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(12000ft ; INaab fc Burkertl (|2003ft ; Uesseit et all (|2005ft and 
ICox et alT i 20061 ) studied the kinematics and photomet- 
ric shapes of simulated merger remnants in order to un- 
derstand the origin of boxy and disky isophotes and the 
connections between the kinematic properties and the 
gas fraction of the merger event, INaab et all (|2006aft 
showed that the line-of-sight vel ocity dispersions are 
stron gly infl uenced by gas ( see a lso iHoffman et ahll2009l 
12010ft . while lQu et ail (|2010ft and lDi Matteo et al.l (|2009ft 
showed that multiple minor mergers slow down the rota- 
tion of an elliptical, while major mergers can speed them 
up. 

Recent resul ts from the Atlas3D survey 

(|Cappellari et al.l 1201 lft . which uses integral field 
spectroscopy to study the properties of early-type galax- 
ies, have revealed that the vast majority of all early-type 
galaxies (82-86 per cent) in their volume-limited local 
galaxy sample is fast rotating, while only very few 
early-typ e galaxies are classifi ed as slow or even non 
rotating (Emscllc m et al.l 12011ft . A significant fraction 
of the slow rotating ellipticals show special features such 
as a counter-rotating core. The non rotating early-types 
are usually found in hig hly overden s e eny ironments 
(|Krainovic et al.l 1201 lft . iBois et al.l (|2011ft showed, 
comparing simulations of isolated merger events to the 
Atlas3D results, that major mergers can reproduce 
all kinds of elliptical galaxies that are fast rotating 
as well as population of flat, slow rotating elliptical 
galaxies with kinematically distinct cores. They fail, 
however, to reproduce very round, massive, slow or even 
non-rotating ellipticals that typically live i n very dense 
environments (see also iBurkert et al1l2008ft . supporting 
the idea that the round slow rotators and the fast 
rotating systems actually are two distinct families of 
early-type galaxies with different formation histories. 

Elliptical galaxies consist mainly of a stellar and a 
dark matter component that might provide interest- 
ing information about their formation history. For ex- 
ample, it is long known from observations that there 
exists a dark halo - disk conspiracy for spiral galax- 
ies, i.e., the rotation curves of spiral galaxi es are flat 
(jEinasto et al.lll974t iFaber fc Gallagherill979ft. Detailed 
dynam ical modeling, for e xamp le by iKronawitter et al.1 
(|2000ft and lGerhard et al.l (|2001ft . has demonstrated that 
a similar conspiracy might also exist between the dark 
halo and the spheroid of massive, slow rotating elliptical 
galaxies. 

Further studies, using planetary nebulae as tracers for 
the outer dark halo, however have revealed a bimodal- 
ity in the kinematic structures of the outer regions of 
elliptical galaxies: Some lower-mass ellipticals show de- 
clining velocity dispersion profiles, similar to a Keplerian 
mass distribution, indicating a constant m ass-to-light ra- 
tio a n d thus a shallow dark ma tter halo (jMendez et al.1 
12001b iRomanowskv et all 12003) , while especially more 
massive, slow rotating galaxies show flat dispersion pro- 
files as expected for a classical, extended da rk matter 
halo (|Napolitano et al.ll2001b iPeng et al.ll2004ft . The ori- 
gin of this bimodality has been discussed, suggesting for 
example that extremely elongated orbits of stars and 
planetary nebulae in the outskirts of low-ma ss ellipticals 
can m ask even a massive dark matter halo (jDekel et al.1 
12005ft . 

A more detailed look at the de-projected dark mat- 



ter component of ellip tical galaxies i nside the ha lf-light 
radius is provided bv iThomas et al.l (120071 12009). using 
Schwarzschild modeling of the observational kinematic 
and photometric data of early-type galaxies in the Coma 
cluster. They show that the best fitting models contain 
10 to 50 per cent dark matter within the half-light radius. 
Furthermore, their models without dark matter halos are 
not able to fit the observations. 

Another approach to study the dark matt er content of 
ellipti cal galaxies comes from strong lensing. lAuger et al.l 
(2010), studying 73 early-type galaxies from the SLACS 
survey, found that the total dark matter plus stellar den- 
sity profiles of these ellipticals can be described on av- 
erage by a power law with a slope of —2.078, with the 
steepness of the individual slopes correlating with the ef- 
fective radius and the central density of t he stellar com- 
ponen t. A similar result was found by Barnabe et al. 
(|2011ft for 16 early-type galaxies from SLACS, where 
they combined the constraints from gravitatio nal lensing 
with s tella r kinematics. Furth er studies by R uff et all 
(|2011ft and lBolton et al.l (|2012ft revealed that this is also 
the case f or strong-lensing e arly-type galaxies at higher 
redshifts: iRuff et al.l (|2011ft found that their 11 early- 
type galaxies in a redshift range of 0.2 < z < 0.65 can 
be fit on average by a slope of —2.16, while iBolton et aTl 
(2012) found for their sample of 79 early-type galaxies 
in the redshift range of 0.1 < z < 0.6 a density slope of 
on average —2.11. Both report a slight trend to flatter 
slopes at higher redshifts. 

An additional method to measure the dark matter con- 
tent of early-type galaxies is weak lensing, which enables 
measurements to mu ch larger radii than strong lensing. 
iGavazzi et al.l (j2007f ) showed for 22 early-type galaxies 
from weak lensing that their total density profiles can be 
described by a power law with a slope of approximately 
—2 for radii as large as 300 kpc, which is up to 100 ef- 
fective radii. 

In a recent study iLvskova et al.l ((2012) investigated 
a sample of cos mological simulati ons of massive galax- 
ies presented by lOser et ail (|2010ft and found a remark- 
able uniformity of the present day isothermal total stel- 
lar mass profiles in good agreement with lensing results. 
In this work we investigate the structure of dark mat- 
ter halos (density distribution and kinematics) around 
simulated ellipticals in more detail and relate it to the 
formation mechanisms and histories of elliptical galax- 
ies. Using a variety of different simulations allows us 
to study the effects of the multiple merger evolution in 
comparison to the major merger scenario as well as the 
influence of the environment on the resulting dynamical 
profiles. 

The paper is structured as follows: In section [5] we de- 
scribe the simulations used for this work. Section [3] ex- 
plains how we selected and analyzed the elliptical galax- 
ies from the simulations. In section |4] we discuss our 
results. 

2. NUMERICAL SIMULATIONS 

We use ellipticals formed in different scenarios: Iso- 
lated binary merger simulations with controlled ini- 
tial conditions (hereafter Binary Merger), hydrodynam- 
ical cosmological zoom-in simulations (hereafter Cosmo- 
Zoom) and large-scale hydrodynamical cosmological sim- 
ulations (hereafter Magneticum Pathfinder). All Simula- 



The Dark Halo - Spheroid Conspiracy in Elliptical Galaxies 



3 



tions were performed using extended versions of the par- 
allel TreePM-SPH-code GADGET-2 (lSpringelll200"5T ) and 
include simplified merger scenarios as well as simulations 
within the cosmological framework, where halos merge in 
a hierarchical fashion. 

GADGET-2 is bas ed on an entropy-conservin g for- 
mulation of SPH (jSpringel fc Hernquistl l2002ft and 
in its standard version includes radiative cooling 
for a primordial mixture of hydrogen and helium 
(|Katz et al.l Il996t ). Star formation and the asso- 
ciated supernova feedback is included usin g a sub- 
resolution model (jSpringel fe Hernquistl 120031 ) and as- 
sume s a Salpeter initial mass function (IMF) (jSalpeterl 
Il955h . In this model the ISM is treated as a 
two-phase medium dMcKee fc Ostrikerlll977t [Efstathiou 
l2000t [Johansson "fc~ Efstathiou 2006) in which cold clouds 
are embedded in a tenuous hot gas at pressure equilib- 
rium. In the cosmological simulations this also includes 
heating by t he time dependent but sp atially uniform UV 
background ([Haardt fe Madaulll996l ). 

The Magneticum Pathfinder simulations also follow the 
pattern of metal pr oduction from the past histo ry of cos- 
mic star formation ([Tornatore et al.ll200l |2007| ) . This is 
done by computing the contributions from both Type II 
and Type la supernovae, and energy feedback and metals 
are released gradually in time, with the appropriate life- 
times of the different stellar populations. This treatment 
also includes, in a self-consistent way, the dependence of 
the gas cooling on the local metallicity. The feedback 
scheme in this case assumes a Chabrier IMF (jChabrierl 
l2003f) . It additionally includes kin etic feedback mimick- 
ing t he effect of supernova winds ([Springel fc Hernquistl 
|2003| ) with its parameters fixed to a wind velocity of 
« 250 kms" 1 . 

Black Holes (BH) are included in the binary merger 
simulations and the Magneticum Pathfinder simulation 
as sink particles. BH feed back is taken into acco unt ac- 
cording to the model from lSpringel et a"H ([2005bl ). where 
the BH sink particle accretes gas from the surround- 
ing medium according to a Bondi-Hoyle accretion model, 
limited to the Eddington limit, and gives back thermal 
feedback to the surrounding medium. Two BHs are as- 
sumed to merge instantly as soon as they enter each oth- 
ers smoothing length and their relative velocity is below 
the local sound speed. 

2.1. Binary Merger Simulations 

The classical formation scenario for an elliptical galaxy 
is the major merger scenario, where we artificially set up 
two spiral galaxies and collide them on specified orbits. 
We analyze a set of 10 high resolution ellipticals formed 
in such a major merger scenario, which we will refer to as 
Binary Ellipticals hereafter. For a detailed description of 
the simulations used in this work, especi ally the details of 
the setup of the progenitor sp i rals, se e Uohansson et al.1 
(|2009af) and Uohansson et al.l (|2009b[ ). Our sample of 
simulations consists of 4 spiral-spiral mergers with a mass 
ratio of 1:1 for the progenitor galaxies, 5 spiral-spiral 
mergers with a mass ratio of 3:1 and one mixed merger. 
For the mixed merger we collide a spiral galaxy with an 
elliptical which is a remnant of a 3:1 spiral-spiral merger 
itself. 

The ga laxies were set up foll owing the method pre- 
sented in ISpringel et al.l ((2005b) . The dark halo virial 



velocity for all primary galaxies except 31OBH2_09_320 
is Wvh- = 160 kms" 1 , halo 31OBH2_09_320 is set up with 
w v i r = 320 kms -1 . Here v v - u - defines the dark matter 
virial mass and virial radius of the spiral galaxy: 



M vir = 



c;: 



IQGHq 

Wvir 



lOi?, 



(I 



(1) 

(2) 



The disks are set up with different initial gas fractions, 
/gas = 0.0, / gas = 0.2 and / gas = 0.8, with the rest be- 
ing disk stars. Simulation 110BHNBO_13 is a simulation 
without a bulge and a gas fraction of / gas = 0.0. All 
of the simulations, except 31ASF2_13 and 31ASF2.13, 
also include Black Holes (BH) as sink particles and BH 
feedback. 

We adopt three orbital geometries, G13, G09 and G01, 
according to iNaab fc Burkertl ((2003). Since we were in- 
terested in studying the effects of the initial gas frac- 
tion and the merger ratios on the merger remnant, we 
choose for the majority of our simulations the same or- 
bital configurations, the G13 orbit. This orbits geometry 
corresponds to an inclination of i\ = —109 and a percen- 
ter argument of oj\ — 60 for the first progenitor galaxy, 
i.e., in case of an unequal mass merger the more massive 
galaxy, and i<i — 180 and ui2 = for the second pro- 
genitor galaxy. Both galaxies approach each others on 
parabolic orbits, with the merger taking place at about 
1.5 Gyr after we started the simulation. The simulations 
were evolved to 3 Gyr. The G09 and G01 orbits were in- 
cluded to understand the impact of the orbit on the rem- 
nant. For the G09 orbit, the parameters are i\ = —109 
and uii = for the first progenitor galaxy and in = 180 
and UJ2 = for the second progenitor galaxy, the G01 or- 
bit matches a geometry of i\ = 0, uj\ = and in = 180, 
uj-2 = 0. All three orbits are parabolic, with a pericentric 
distance of r pcr i = r^.i + r^2, with r^i and r^,2 the disk 
scale radii for the first respective the second progenitor 
galaxy. G13 and G09 are orbits where the orientations of 
the progenitors disks are not in the orbital plane, while 
G01 is a retrograde orbit with the progenitors being in 
the orbital plane . For more details on these orbits see 
INaab fc Burkertl ((2003). Table [TJ contains a summary of 
all simulation parameters used for this study. 

The mass resolution for all binary merger simulations 
is Mdm = 2.25 x 10 6 A/ Q ft- _1 for dark matter particles 
and M gas = M sta rs = 1-30 x 10 5 M Q h~ 1 for gas and star 
particles. The gravitational softening length was set to 
edm = 0.083/i~ 1 kpc for dark matter particles and to 
e gas = £stars = 0.02/i _1 kpc for gas and star particles. 
For all simulations we used h = 0.71 and a baryonic 
mass fraction of f^B = 0.044. 

2.2. Cosmological Zoom-In Simulations 

To study the formation of elliptical galaxies within the 
cosmological context, i.e., from multiple mergers, we an- 
alyzed a set of galaxies extracted from cosmological sim- 
ulations of structure formation. Our sample consists of 
17 zoom-in re-simulations of individual halos, hosting a 
central spheroidal galaxy at z=0. Additionally, we in- 
cluded four companion ellipticals, i.e., massive spheroids 
that are substructures within larger halos. We will refer 
to the central spheroidal galaxies as CosmoZoom Ellip- 
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TABLE 1 

Binary merger simulation sample at a timestep of 3 Gyr 



Model 


Ratic-W 


Orbit( b ) 


f (c) 


M Gal W 


M DM (C) 




jncw(g) 


Rl/2 W 


fDM (i) 


f0.5 (k) 
J DM 


11 OBH2 13 


1:1 


G13 


0.2 


1.84 


4.56 


712 870 


10.71 


4.67 


0.24 


0.14 


11 OBH2 09 


1:1 


G09 


0.2 


1.84 


4.53 


713 997 


11.21 


4.45 


0.23 


0.13 


11 OBH0 13 


1:1 


G13 


0.0 


1.92 


4.39 


742 490 


0.00 


6.14 


0.29 


0.18 


11 OBHNB0 13 


1:1 


G13 


0.0 


1.42 


4.49 


548 320 


0.00 


7.56 


0.46 


0.32 


mix 11 OBH2 13 


1:1 


G13 


0.2 


2.09 


4.77 


809 734 


11.79 


5.34 


0.26 


0.15 


31 OBH2 13 


3:1 


G13 


0.2 


1.20 


3.22 


466 326 


9.34 


4.75 


0.29 


0.17 


31 OBH8 13 


3:1 


G13 


0.8 


1.07 


3.31 


416 201 


50.84 


2.40 


0.15 


0.08 


31 ASF2 13 


3:1 


G13 


0.2 


1.21 


3.22 


469 070 


9.77 


4.53 


0.28 


0.16 


31 ASF2 01 


3:1 


G01 


0.2 


1.23 


3.24 


474 622 


11.63 


3.96 


0.24 


0.13 


31 OBH2 09 320 


3:1 


G09 


0.2 


9.62 


25.60 


465 699 


10.14 


9.37 


0.28 


0.17 



disks of the progenitor galaxies; (d) Stellar masses within 10 per cent of the dark matter virial radius in 10 11 AfQ; (e) Dark matter masses 
within 10 per cent of the dark matter virial radius in 10 11 AfQ; (f) Number of stellar particles within 10 per cent of the dark matter virial 
radius; (g) Fraction of newly formed stars since the start of the simulation in comparison to the total number of stars at the final output 
time t = 3 Gyr; (h) Effective radius of the stellar component of the galaxy, calculated as three dimensional half-mass radius, in kpc (i) 
Fraction of dark matter relative to the stellar component within the half-mass radius (k) Fraction of dark matter relative to the stellar 
component within 0.5-Ri/2 

TABLE 2 

CosmoZoom Ellipticals sample at z=0 



Model 




M Gal t b > 






yncw (c) 




f («) 

J in— situ 


fDM W 


r0.5 (i) 

J DM 


0040.2 


36.08 


6.94 


26.42 


84 786 


2.02 


12.91 


23.07 


0.43 


0.29 


0053.2 


22.58 


9.65 


26.73 


117 809 


1.43 


13.03 




0.37 


0.23 


0069.2 


25.12 


6.86 


18.44 


83 804 


1.49 


8.84 


21.78 


0.31 


0.20 


0089.2 


14.94 


7.27 


15.78 


88 772 


0.57 


10.43 


16.43 


0.34 


0.22 


0094.2 


13.46 


6.65 


17.36 


81 258 


0.80 


7.53 


25.99 


0.27 


0.15 


0125.2 


12.03 


6.02 


15.36 


73 546 


0.90 


9.08 


22.45 


0.33 


0.20 


0162.2 


8.72 


5.06 


11.70 


61 808 


2.11 


9.78 


12.92 


0.40 


0.27 


0163.2 


9.78 


4.89 


9.96 


59 701 


0.90 


10.38 


14.96 


0.37 


0.27 


0175.2 


9.59 


5.11 


13.80 


62 401 


1.05 


7.37 


26.96 


0.31 


0.20 


0190.2 


8.09 


4.37 


8.11 


53 401 


2.72 


6.99 


14.67 


0.29 


0.20 


0204.2 


8.15 


3.73 


8.27 


45 554 


1.63 


6.50 


15.63 


0.25 


0.15 


0204.4 


7.96 


4.07 


7.48 


398 025 


1.40 


6.36 


23.86 


0.20 


0.11 


0215.4 


8.37 


4.36 


9.30 


425 565 


1.50 


4.25 


36.14 


0.19 


0.10 


0408.4 


3.96 


2.30 


3.19 


224 648 


0.94 


4.51 


20.88 


0.16 


0.10 


0501.4 


3.77 


2.33 


4.42 


227 403 


2.70 


3.27 


39.31 


0.19 


0.10 


0616.4 


3.50 


2.38 


3.72 


232 935 


4.00 


5.02 


29.90 


0.27 


0.17 


0664.4 


2.93 


1.70 


3.30 


166 379 


2.46 


2.81 


38.12 


0.16 


0.08 



Note. — (a) Total stellar mass within the dark matter virial radius in 10 Mq; (b) Stellar mass within 10 per cent of the dark matter 
virial radius in 10 11 Mq; (c) Dark Matter mass within 10 per cent of the dark matter virial radius in 10 11 Mq; (d) Number of stellar particles 
within 10 per cent of the dark matter virial radius; (e) Fraction of newly formed stars since a redshift of 0.27, which is approximately 3 Gyr, 
for comparison with the major merger sample; (f) Effective radius of t he stellar compon ent of the galaxy, calculated as three dimensional 
half-mass radius, in kpc (g) Fraction of stars formed in-situ taken from Oser ct al. (2012) in per cent, (h) Fraction of dark matter relative 
to the stellar component within the half-mass radius (i) Fraction of dark matter relative to the stellar component within 0.5i?!/2 



TABLE 3 

CosmoZoom companion sample at z=0 



Model 


AW aJ 


A/DM lb) 


N G ^) 


jnew (d) 


Ri/2 (a) 




/DM (g) 


f 0.S (h) 
J DM 


0040.sl 


4.55 


14.07 


55 581 


2.05 


6.62 


445.0 


0.28 


0.18 


0069_sl 


1.84 


4.96 


22 451 


1.15 


4.31 


317.5 


0.30 


0.20 


0069.s2 


2.32 


7.42 


28 393 


1.52 


4.54 


339.5 


0.30 


0.20 


0069.s3 


1.97 


5.16 


24 013 


3.72 


3.84 


340.0 


0.29 


0.19 



Note. — (a) Stellar mass within 10 per cent of the dark matter virial radius in 10 11 A/q; (b) Dark Matter mass within 10 per cent of 
the dark matter virial radius in 10 11 Mq; (c) Number of stellar particles within 10 per cent of the dark matter virial radius; (d) Fraction 
of newly formed stars since a redshift of 0.27, which is approximately 3 Gyr, for comparison with the major merger sample; (e) Effective 
radius of the stellar component of the galaxy, calculated as three dimensional half-mass radius, in kpc; (f) Maximum circular velocity at 
infall into the parent halo (g) Fraction of dark matter relative to the stellar component within the half-mass radius (h) Fraction of dark 
matter relative to the stellar component within 0.5i?i/2 
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TABLE 4 
Magneticum BCGs at z=0 



Model 


J^tot(a) 


AW b) 






Rl/2 (C) 


fDM (l) 


r().5 (g) 

J DM 


00 


174.4 


52.2 


508.6 


142 140 


62.68 


0.71 


0.52 


02 


178.6 


69.4 


439.8 


163 955 


40.47 


0.51 


0.26 


03 


152.8 


88.3 


359.4 


182 988 


30.98 


0.30 


0.14 



Note. — (a) Total stellar mass within the dark matter virial radius in 10 11 Mq\ (b) Stellar mass within 10 per cent of the dark matter 
virial radius in 10 11 Mq; (c) Dark Matter mass within 10 per cent of the dark matter virial radius in lO n M0; (d) Number of stellar 
particles within 10 per cent of the dark matter virial radius; (e) Effective radius of the stellar component of the galaxy, calculated as three 
dimensional half-mass radius, in kpc. (f) Fraction of dark matter relative to the stellar component within the half-mass radius (g) Fraction 
of dark matter relative to the stellar component within 0.5i?i/2 



ticals hereafter, and to the substructures as CosmoZoom 
Companions. 

The parent cosmological box of 72 3 /i -3 Mpc 3 was sim- 
ulated using 512 3 DM particles with a particle mass of 
Mdm = 2x10 8 Mq/i -1 and a comoving gra vitational soft- 
ening length of 2.52ft" 1 kpc. A WMAP3 (|Spergel et al.l 
2007), ACDM cosmolo gy was adopted, with erg = 0.77, 
Oa = 0.74, il m — 0.26 and h = 0.72 and an initial slope 
for the power spectrum of n s — 0.95. 

From this dark matter only simulation, halos of dif- 
ferent masses, ranging from 10 11 M Q /i _1 to 10 13 Mg/i -1 , 
were selected at z — 0. All dark matter particles closer 
than 2 x T2oo to the halo center at any snapshot are 
traced back in time. We replace all dark matter particles 
identified that way by dark matter and gas particles at 
higher resolution with J7dm = 0.216 and JIb = 0.044. 
The details of the re-simulation method are described in 
lOser et al.l IpOinh . 

These simulations use cooling for a primordial gas com- 
position and star-formation but do not include any black 
hole treatment. The initial co nditions were crea ted us- 
ing GRAFIC and LINGERS (|Bertschingerl l200l . The 
simulations were evolved from z ~ 43 to z = 0. 

To achieve a proper resolution even for the smaller 
halos we performed the re-simulations at two different 
resolutions. The most massive halos were re-simulated 
with twice the spatial resolution of the original dark 
matter only box. In these re-simulations the parti- 
cle masses are Mdm = 2.1 x 10'Mg/i -1 and Af gas = 
Mtars = 4.2 x lOpM^hr 1 with gravitational softening 
set to £dm = 0.89/i _1 kpc and e gas = e s t a rs = 0Ah~ kpc 
for dark matter, gas and star particles, respectively. To 
study the effects of the gas physics and the stellar compo- 
nent on the dark matter, the halos of this re-simulation 
level were also re-simulated with dark matter only at the 
same resolution. 

For the less massive halos we used four times the spa- 
tial resolution of the original box, and particle masses 
of Mdm = 3.6 x 10 6 M Q /i _1 for the dark matter par- 
ticles and M gas = M stara = 7.4 x 10 5 M Q /i _1 for the 
gas and star particles with gravitational softening set to 
6dm = 0.45/i _1 kpc and e gas = e stars = 0.2/i _1 kpc re- 
spectively. Table [3] contains all ellipticals extracted from 
the re-simulations used for this study, and Table [3] con- 
tains the companion halos, labeled with the name of their 
host central halo followed by an S and the number of the 
substructure. 

2.3. Magneticum Pathfinder Simulation 

We also extracted central galaxies of clusters from a 
hydrodynamical, cosmological simulation (Magneticum 



Pathfinder, Box3/hr, Dolag et al 2012, in preparation). 
These galaxies were extracted from a 128 3 /i -3 Mpc 3 
Box simul ated using 2 x 576 3 particles and adapting a 
WMAP7 (jKomatsu et al.l 1201 ID ACDM cosmology with 
a 8 = 0.809, h = 0.704, fi A = 0.728, n M = 0.272 and 
Ob = 0.0456 and an initial slope for the power spec- 
trum of n s = 0.963. Dark matter particles have a mass 
of Mdm = 6.9 x 10 8 M Q /i _1 , gas particles have a mass 
of approximately Mcas = 1.4 x 10 8 M Q /i~ 1 , depending 
on their enrichment history, and stellar particles have 
about Mg tars = 3.5 x 10 7 Mq/i _1 , depending on the state 
of the underlying stellar population. For all components 
the gravitational softening length is set to 3.75/i _1 kpc. 
In this simulation, one gas particle can form up to four 
stellar particles. 

This simulation follows the gas using a l ow-viscosity 
SPH scheme to properly track turbulence (jDolag et al.1 
|2005| ) and the cooling and star-formatio n description fol- 
lows the pa ttern of metal production ([Tornatore et al.l 
2004 I2007T) . including in a self-consistent way the de- 
pen dence of the gas co oling on the local metallic- 
ity (jWiersma et al.l 120091 ) . It al so includes the feed- 
back caused by galactic winds (|Springel fc Hernquistj 
2003T) as well as feedb ack from black holes according to 
Soringel et all (|2005bh . 

We selected the three halos among the most massive 
ones of which the central galaxy contained at least 10 5 
star particles and was not currently undergoing a merger. 
We will refer to these massive central cluster galaxies as 
Magneticum BCGs hereafter. Table [4] summarizes the 
properties of all three ellipticals. 

3. RESULTS 

We identify all stars within 0.1i? v j r to belong to the 
central elliptical, similarly to lOser et ail (|2010l ). inde- 
pendent of the simulation type. The galaxy masses used 
throughout this work are calculated according to these 
definitions. 

For the companion ellipticals from the cosmological 
simulations, where significant parts of the dark matter 
have been stripped already and thus the virial radius 
cannot be obtained from the dark matter mass of the 
subhalo, we use the maximum circular velocity v max to 
calculate the virial mass and the virial radius the halo 
had before the infall into the parent halo, using equa- 
tions [1] and [2] with u max as u v ir- These companion ellip- 
ticals are particularly interesting since our observational 
compar ison sample, the elli ptical galaxies in the Coma 
cluster (jThomas et alJl2C)07f ) , are embedded in the larger 
Coma cluster environment themselves and thus suppos- 
edly have a similar formation scenario to the companion 
ellipticals. 
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For most of this work we use the elliptical galaxies at 
the final output time of 3 Gyr for the binary merger simu- 
lations and the ellipticals identified at a redshift of z = 
for all cosmological simulations. First we calculate the 
three dimensional stellar half-mass radius i?i/ 2 for all el- 
liptical galaxies in our sample, i.e., the radius which con- 
tains 50% of all stars that were identified to belong to the 
galaxy. Next the intrinsic velocity dispersion and density 
within spherical equal-mass shells around the center of 
our simulated galaxies are calculated for stars and dark 
matter separately as well as for both components com- 
bined: Within each shell we first calculate the intrinsic 
velocity dispersion for stars and dark matter separately, 
as 



(3) 



with oy = y/ (v 2 ) — (v r ) 2 and erg and 0$ analogous. Since 
dark matter and stellar particles have different masses 
within all our simulations, we need to calculate the total 
intrinsic mass-weighted velocity dispersion for the com- 
bined profiles within each shell as 



Ctot 



' m DM o-p M + msug 
m D M + m 3 



(4) 



to obtain a total velocity dispersion that is representa- 
tive for the whole potential. Profiles for three example 
halos are shown in Figure [1] for the density in the upper 
panels and the velocity dispersion in the lower panels. 
We choose one example halo for each simulation type, 
for the Binary (left), the CosmoZoom (central) and the 
Magneticum (right) Ellipticals. 

A power law is fit to these velocity dispersion and den- 
sity profiles, to stars and dark matter separately as well 
as to the combined profiles. The fit to the combined 
profiles is shown as a green dashed line in Figure [1] for 
the three example halos. We were especially interested 
in the transition area, where the dark matter part be- 
comes dominant in its contribution to the total density 
and velocity profiles. This transition area is character- 
ized well by the stellar half- mass radius i?i/ 2 for all galax- 
ies, and thus we choose the fitting range to depend on 
Ri/2- For the Binary Ellipticals and the CosmoZoom El- 
lipticals the best-fitting power law index for the velocity 
dispersions is determined in the radius regime 0.3i?iy 2 
to 2i? 1 / 2 , the power law for the density is fit between 
0.3i?i/2 and 4i?!/ 2 . For the CosmoZoom Companions 
and the Magneticum BCGs the best-fitting power laws 
are determined with a lower limit of 0.5i?!/ 2 instead of 
0.3i?i/2 due to resolution limits. The lower limits are 
chosen to ensure that the innermost particles included 
are at least 3 times the smoothing length away from the 
center. The upper limits are chosen due to comparability 
with o bservations of Coma ellipticals by iThomas et al.1 
(|2007f ). where density slopes are available up to 4i? e ff 
while velocity dispersion slopes are only available up to 
2-Roff- 

Figure [2] shows the combined dark matter and baryonic 
density profiles multiplied by r 2 for all Binary and Cos- 
moZoom Ellipticals normalized to the density at 0.3i?i/ 2 
and for all CosmoZoom Companions and Magneticum 
BCGs normalized at 0.5i?!/ 2 , illustrating that a power 
law fit is actually a good approximation for most of our 



ellipticals. The only exceptions are the halos with an 
unusually dominant central stellar component, which are 
the least massive of the CosmoZoom Ellipticals with four 
times the spatial resolution and some of the CosmoZoom 
Companions. These halos have slightly curved combined 
density profiles and the best fitting power laws have gen- 
erally steeper slopes. 

Figure [T] and Figure [2] also show that the power law fit 
to the total density profiles of the CosmoZoom Ellipticals 
and the Magneticum BCGs is a good approximation to 
even larger radii than 4i? 1/ / 2 , even though the profiles be- 
come less smooth due to the presence of satellite galaxies. 
This is in agreement with o bservations from weak lens- 
ing by iGavazzi et al.l ([2007) , who found for 22 central 
elliptical galaxies that their total density profiles can be 
fit on average by a power law with a slope of 7 ~ —2 
over two decades in radius. A similar behavior can also 
be seen for the CosmoZoom Companions, although their 
density profiles break down at about 10i?i/ 2 due to the 
fact that their dark matter halos got stripped during the 
infall into the larger parent halo. For the Binary ellip- 
ticals, the density profiles already break down at about 
5i?i/ 2 due to the fact that the simulations do not reach 
further out. 

3.1. Density and velocity dispersion slopes 

For a spherical isothermal system the solution of the 
Jeans Equation is 



p(r) 



2nGr 2 ' 



(5) 



More generally, under the assumption that both the ve- 
locity dispersion and the density can be described by 
simple power laws, i.e., p(r) = Ar 1 and <r(r) = Br? with 
A and B being constants, the Jeans Equation has the 
following solution: 



p(r) 



Ca(r) 



(6) 



AirGr 2 

with C being a dimensionless constant. Thus, the slopes 
of the dispersion and the density are correlated by 



p = O.57 + I. 



(7) 



This correlation is shown in both panels of Figure [3] and 
Figure S] as a solid black line and in Figure [5] and Figure [S] 
as dotted black line. These relations hold for all systems 
with constant anisotropy as a function of radius, since 
different values for the constant anisotropy only change 
the constant C and thus Equation [7] does not change. 

Figure [3] shows the total density slopes and total ve- 
locity dispersion slopes for all Binary Ellipticals in the 
upper panel and for all CosmoZoom Ellipticals and Mag- 
neticum BCGs in the lower panel. We can clearly see that 
all total slopes of the Binary Ellipticals lie close to the so- 
lution of the Jeans Equation, i.e., they have total density 
slopes around 7tot — —2.1 and total velocity dispersion 
slopes around /3 to t = 0, implying that these ellipticals are 
fairly close to spherical systems with constant anisotropy, 
or even isotropic in some cases. 

For our sample, neither the choice of different orbits 
(G01, G09 or G13) nor the merger type (i.e., if it is a 
3:1, a 1:1 or an E-SP merger) or the presence of a Black 
Hole changes the resulting density slopes significantly, 
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Fig. 1. — Density (upper panel) and velocity dispersion (lower panel) profiles for three example halos. Left panels: 3:1 Binary Elliptical 
31_OBH2_13. Central panel: CosmoZoom Elliptical 190. Right panel: Magneticum BCG 03. Red solid line: stellar profile; blue solid line: 
dark matter profile; black line: combined profile; green line: power law fit to the combined profile. Dashed black lines for the left and 
central panel: 0.3i?!/2 an d &R\/2 f° r the density, 0.3i?i/2 an d 2Ri/2 f° r the velocity; for the right panel:0.5Ri/2 an d ^R\I2 f° r the density, 
0.5-R]y2 an d 2Ri/2 f° r the velocity. The arrows mark the position of the half-mass radius. The dash-dotted blue line in the central panel 
shows the dark matter profile for the dark matter-only re-simulation of CosmoZoom Elliptical 190. 
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Fig. 2. — Combined dark matter and baryonic density profiles for all ellipticals, multiplied by r 2 and normalized to the density value 
at 0.3i?i/2 (0-5-R]/2 for the CosmoZoom Companions and the Magneticum BCGs). From left to right: Binary Ellipticals, CosmoZoom 
Ellipticals with twice the spatial resolution, CosmoZoom Ellipticals with four times the spatial resolution, CosmoZoom Companions, 
Magneticum BCGs. The dashed lines mark 0.3i? 1; /2 ancl ^Ri/2 f° r the first three panels, 0.5Ri/2 and iRi/2 in the last two panels. 



while the variation of the initial gas fraction causes the 
only significant shift in the density slope. The simulation 
with 80% initial gas fraction has a significantly steeper 
total density slope than the comparable simulation with 
20% initial gas fraction. At present day, mergers with 
high gas fractions are unlikely, but at a higher redshift 
gas rich major mergers are much more frequent. Still, our 
80% gas merger has very large and extended progenitor 
gas disks, which is unrealistic for high-z disks. 
The CosmoZoom Ellipticals show a much larger vari- 



ety of total density slopes (see the lower panel of Figure 
[3]) than the Binary Ellipticals, although they also have 
flat total velocity dispersion curves, i.e., the slopes of the 
power law fits are close to zero. This is in agre e ment with 
a detailed analysis presented in[Lyskova ct al. (2012]), us- 
ing the same simulation set as lOser et al.l ()2010l ). While 
the CosmoZoom Ellipticals with the flatter density slopes 
around 7tot = —1-9 and total velocity dispersion slopes 
around /3tot = 0.05 are close to the solution of the Jeans 
Equation and thus fairly close to spherical systems with 
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Fig. 3. — Slopes of the total velocity dispersions /3tot against the 
slopes of the total density profiles 7 tot . Upper panel: Results from 
the Binary Ellipticals: 1:1 spiral merger (black), 3:1 spiral merger 
(blue), Elliptical-spiral merger (cyan) and 3:1 spiral merger with 
80% gas (pink). Lower panel: Results from the cosmological sim- 
ulations: CosmoZoom 2X Ellipticals (bright green), CosmoZoom 
4X Ellipticals (dark green), CosmoZoom Companions (yellow) and 
Magneticum BCGs (violet). For all ellipticals the errors are RMS- 
deviations to the fit. Black line: analytic solution for a spherical 
system with constant anisotropy. 

constant anisotropy, the CosmoZoom Ellipticals with the 
steeper density profiles are not. The steeper the total 
density slopes are, the larger is the deviation from the 
Jeans solution. This could be due to a combination of 
gradients in the anisotropy and non-spherical effects that 
are not included in the simple spherically symmetric ap- 
proach of this paper. The details will be discussed in 
a subsequent paper (Remus et al, in preparation). We 
also find that there is no difference between the behav- 
ior of the CosmoZoom Ellipticals and the CosmoZoom 
Companions. All three Magneticum BCGs are close to 
the Jeans solution and have relatively fiat total density 
slopes compared to the majority of the CosmoZoom El- 
lipticals, with their total velocity dispersion slopes are 
flat as well. 

Figure [4] shows the stellar velocity dispersion slopes 



Fig. 4. — Same as Figure \3\ but for stellar velocity disper- 
sion slopes versus total density slopes. Red triangles show the 
total density and stellar velocit y dispersion slopes fo r the Coma 
early-type galaxies presented in Thomas ct al. (2007). The yel- 
low (upper panel) respective blue (lower panel) triangle shows the 
slopes for the ma s sive s trong lensing early-type galaxy studied by 
ISonnen feld ct al. (2012). Dashed lines: linear fits to the Coma 
values (red) and all cosmological simulations (green). 



/3* and the total density slopes 7tot for the Binary El- 
lipticals in the upper and the CosmoZoom and Mag- 
neticum Ellipticals in the lower panel. This figure also 
includes the stellar velocity dispersion and total density 
slopes obtained fro m observations of the Coma ellipti- 
cals as presented in I Thomas et al.l (|2007) and the mas- 
si ve, strong-lensing e arly-ty pe galaxy discussed in detail 
in ISonnenfeld etahl (j2012n . As can clearly be seen in 
the upper panel, our sample of Binary Ellipticals cannot 
reproduce the range of slopes that is seen for the Coma 
cluster ellipticals, and the observations show no cluster- 
ing around the values of the Binary Ellipticals. Thus the 
scenario of a present-day major merger seems unlikely to 
be the dominant formation scenario for the Coma ellip- 
ticals. Binary merger between two high-redshift spirals 
might produce different results, because the initial con- 
ditions for those galaxies would look very different, i.e., 
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Fig. 5. — Same as Figure \3\ but for both velocity and den- 
sity slopes of the dark matter alone. Red circles: slopes of the 
velocity and density from the CosmoZoom dark matter-only re- 
simulations. Dotted black lines: analytic solution for a spherical 
system with constant anisotropy. The density slopes cluster in a 
range of —1.8 < 7dm J; —1.4 excluding the CosmoZoom Compan- 
ions (yellow circles) and the dark matter-only re-simulations, as 
explained in the text. 

they would for example have no large stable gas disks 
and different dark matter halos, thus we cannot exclude 
high redshift binary mergers to be the dominant forma- 
tion mechanism for the Coma ellipticals. 

As can be seen in the lower panel of FigureHJ the range 
of values of the slopes of the CosmoZoom Ellipticals is 
similar to the range found for t he observed Coma cluster 
ellipticals (| Thomas et al . 1120071 ) . although there is a slight 
offset in the velocity dispersion slopes with respect to 
the observations. The Magneticum BCG are in good 
agreement with the slopes found for the more massive 
Coma Ellipticals, although the observational sample does 
not include the BCGs since the data of the Coma BCGs 
reach out to only 0.5i?!/ 2 - 

Interestingly, if we look at the density slopes against 
the velocity slopes for the dark matter component only, 
the values cluster around 7dm = — 1-5 for the Binary El- 
lipticals and around 7dm = —1.67 for the CosmoZoom 



Fig. 6. — Same as Figure |31 but for both velocity and density 
slopes of the stellar component alone. Dotted black lines: ana- 
lytic solution for a spherical system with constant anisotropy. The 
slopes show no clustering in the density, they are spread over a 
range of -3.25 < 7* < -2.4. 

Ellipticals and the Magneticum BCGs, as shown in Fig- 
ure [5] This is in agreement with the o bservational results 
presented by ISonnenfeld et ail ((2012), who found a dark 
matter density slope of 7dm = —1.7 ± 0.2 for their mas- 
sive strong-lensing early-type galaxy. Only the Cosmo- 
Zoom Companions show different values, which is most 
likely due to the fact that the outer parts of their dark 
matter halos have been stripped significantly during the 
infall in the parent halo, causing a steeper slope. 

We also include in the lower panel of Figure[5]the slopes 
from the fits to the profiles of the dark matter only re- 
simulations of the CosmoZoom-2x Ellipticals, to study 
the influence of the stellar component on the dark mat- 
ter halos directly. As can be seen, the density slopes of 
the dark matter only re-simulations are slightly flatter 
(around 7dm only = —1-46) than for the dark halos that 
contain a stellar component, and both their density and 
velocity dispersion slopes are closer to the theoretical so- 
lution for an isothermal sphere. This can also be seen in 
the central panels of Figure [T] which shows the density 
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and velocity dispersion curves of the dark matter compo- 
nent for an example CosmoZoom Elliptical as solid blue 
line and the corresponding curves for the dark matter 
only simulation of the same halo as the blue dashed line. 
These values for the slopes of dark matter halos from 
dark matter only simulations are in agreement with re- 
sults for the slopes of the central density of dark matter 
halos from hi gh resolution dark m atter only simulations 
presented by iMoore et al.l (|1999D . who found slopes of 
7dm only ~ —1-5, which is the same as what w ould be 
expected for NFW-prohles (Na varro et al.llT996[ ) in this 
radius range. 

We see that the presence of the stellar component sig- 
nificantly alters the distribution of the dark matter. If 
the baryonic component is included in the simulation, 
the dark matter halo is denser in the center and thus 
the dark matter density slope is steeper, although for 
both the simulation with and without baryons the den- 
sity converges to the same values at large radii beyond 
approximately 5i?i/ 2 - This is in agreement with results 
presented recently by Uohansson et al.l (|2012f ) and is due 
to the well know effect of adiabatic contraction, i.e., 
the dark matter particles are pulled inwards due to the 
condensation of the gas in the center of t he halo (e.g., 
Jesseit et aLll2002t [Blumenthal et alJI 19861: iGnedin et al.1 
2004 LI201U see however iDutton et al.ll2007l and, regard- 
ing the effects of expansion due to sudden out flows driven 
by supernovae, iPontzen fc Governatol 120121 ) . An even 
stronger effect caused by the presence of baryons can 
be seen for the velocity dispersion profile in the lower 
panel of Figure [TJ While the velocity dispersion of the 
dark matter component of the simulation with baryons 
slightly decreases with larger radii, the dark matter-only 
simulation shows a velocity dispersion that is strongly 
increasing with larger radii up to approximately 2i? 1 / 2 
and thus the power law fit has a positive slope. 

Figure [6] shows the density and velocity dispersion 
slopes for the stellar component, with the Binary Ellipti- 
cals in the upper and the CosmoZoom and Magneticum 
BCGs in the lower panel. For all ellipticals, the den- 
sity slopes of the stellar component are generally steeper 
than the slopes of the dark matter component. While 
the slopes of the stellar component for the Binary Ellip- 
ticals are all around 7, = —2.9, the slopes of the stellar 
component of the CosmoZoom Ellipticals and the Mag- 
neticum BCGS show a different behavior: In contrast 
to the dark matter component, which has a small range 
of density slopes from —1.8 < 7dm < —1.4 excluding 
the Companion ellipticals, the stellar component covers 
a larger range of slopes (—3.25 < 7, < —2.4) and there 
is no correlation between the stellar slope and the type 
of simulation. This implies that the stellar component is 
responsible for the steepness of the total density slope, 
namely a more dominant stellar component in the center 
of a galaxy leads to a steeper total density slope. 

3.2. The influence of gas and star formation on the 
density slope 

Unlike the stellar and dark matter particles, which are 
collisionless, the gas particles can dissipate their energy 
and thus condense in the center of the galaxy. Thus, new 
stars are mostly formed in the central area of the galaxy. 
The more gas that is present in a merger event, the more 
dominant is this newly formed stellar component in the 
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Fig. 7. — Same as Figure[T]for the 3:1 Binary Elliptical with 20% 
initial gas fraction (left panel) and the 3:1 Binary Elliptical with 
80% initial gas fraction (right panel). 



center. 

This effect is nicely demonstrated in the comparison 
between the two Binary Ellipticals which have identical 
initial conditions apart from the initial gas fraction of 
20% respective 80%. Figure [7] shows the density and 
velocity dispersion profiles for both ellipticals. In case of 
the merger with 80% gas fraction we see two effects that 
change the stellar density profile compared to the case 
of the merger with a gas fraction of 20%: First, we see 
a that the overall density is generally higher since more 
stars have been formed over the whole radius, i.e., the 
normalization A in the profile p{r) = Ar 1 is larger, while 
the dark matter profile for both halos did not change 
significantly. Second, the central part of the elliptical 
is much more compact, as a large amount of stars has 
been formed there. Both effects together cause the total 
density profile to be much steeper in case of the 80% gas 
merger. The half-mass radius of the 80% gas merger is 
with i?i/2 = 2.4 kpc just half as large as the half-mass 
radius of the 20% gas merger, while the fraction of stars 
that are formed during the merger is with /" ow = 51% 
much higher than in the 20% gas merger (/" ew = 9%), 
as can be seen in Table [1] 

This correlation between the steepness of the total den- 
sity slope and the half-mass radius can also be seen for 
the CosmoZoom Ellipticals and Companions, as shown 
in the left panel of Figure HI The steeper the total den- 
sity slope of an elliptical, the smaller the half-mass ra- 
dius. We also find a (weaker) correlation between the 
stellar mass of an elliptical galaxy and its total density 
slope 7tot, as shown in the right panel of Figured! and a 
correlation between the total density slope and the dark 
matter fraction within the half-mass radius, as shown in 
Figure [9l As expected, the total density slope is steeper 
the more dominant the stellar component is compared to 
the dark matter in the inner part of the galaxy, i.e., the 
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smaller the fraction of dark matter within the half-mass 
radius. 

In these figures we included t h e res ults for the 
Coma ellipticals by iThomas et al.1 ()2007j ) as well as 
the result s from the SLACS strong lensing s urvey pre- 
sented by I Auger et al.l (12010D . iBarnabe etall (|2011h and 
ISonnenfeld et al.l (|2012t ). and we see that our result are 
in good agreement with the observations. The only ex- 
ceptions are the Magneticum BCGs for which we have no 
observational counterparts, neither in mass nor in half- 
mass radius. 

Also shown in the left panel of Figure [5] are the total 
density slopes against the half-mass radii for the Cosmo- 
Zoom dark matter-only simulations. As seen before, the 
dark matter-only simulations show much flatter density 
slopes than the simulations including baryon physics, i.e., 
adding the baryons steepens the total density profile of 
the halos. We find that the observations clearly favor the 
slopes given by simulations with baryon physics. 

There are a few early-type galaxies in the Coma obser- 
vational sample that have a slope around 7tot ~ — 2 and a 
very low dark matter fraction (see the upper panel of Fig- 
ure [9]). This kind of ellipticals with very low dark matter 
fractions are also p resent in the strong lensing sample of 
lAuger et al.l ()2010l ) (see the lower panel of Figure ED, bu t 
not in the strong lensing sample of lBarnabe et al.l (2011). 
We cannot reproduce these early-type galaxies with any 
of our simulated scenarios, not even with the Cosmo- 
Zoom Companions. Those early-types seem to have a 
dominant stellar component, but a relatively flat density 
slope. 

In case of the Com a Cluster ellipticals we know from 
IThomas et al.l (|2011l) that these ellipticals with low dark 



matter fractions have large dynamical mass-to-light ra- 
tios in the case of the Coma Cluster ellipticals compared 
to a Kroupa IMF. This means that these ellipticals ei- 
ther have a bottom-heavy stellar initial mass function or 
that their dark matter density is nearly identical to the 
density of the stellar component. 

There have been several recent papers indicating es- 
pecially in case of massive early- typ e galaxies, that 
the IMF is not universal but variable (|Cappellari et all 
2012tlFerreras et alll20ll Ivan Dokkum fc Conrovll2011l: 
Conrov fc van Dokkuml 120121 : Ivan Dokkum fe Conroyl 
2012tlTreu et al.ll2010() . This is interesting since the pre- 
dicted dark matter fractions strongly depend on the as- 
sumed IMF. For example, in case of an IMF like Kroupa, 
the observed dark matter fractions for the Coma Clus- 
ter ellipticals would be much higher (between 40-70%, 
see IThomas et all 1201 11 1 and thus fit quite well to our 
results from the simulations . H owever, as shown b y 
IConrov fc van Dokkuml (|2012l ) and lWegner et all(|2012l) . 
not even a variable stellar IMF can always explain the 
high dynamical mass-to-light ratios that are observed. 

On the other hand, if the dark matter density follows 
the stellar component closely enough, both components 
become indistinguishable from each other and thus the 
stellar mass becomes overestimated in the Schwarzschild 
modeling. This is explained in de tail for the early-typ e 
galaxies in the Abell 262 cluster in lWegner et al.l (f2012j) . 
One possible way to explain an increase of the dark 
matter is by adiabatic contraction, as discussed before- 
hand. Nevertheless, the contraction would have to be 
very strong, stronger then what is seen in our simula- 
tions. Another process that could cause similar stellar 
and dark matter densities is violent relaxation, which 
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has been discussed in lWegner et al.l ([2012) as well. 

It is also possible that our simulations simply do not in- 
clude ellipticals that are the equivalents to those observed 
early-types. It is possible that these early-type galaxies 
are actually spirals that suffered from tidal and gas strip- 
ping while they entered a dense environment lik e, for ex- 
ample, the Coma cluster. iGunn fc Gottl (fl972T) already 
suggested ram pressure stripping to be an efficient way 
to form SO early-type galaxies in dense environments, 
and ongoing strip ping has been recently observed for the 
Virgo Cluster by ([Abramson et al.ll2011l) , but we have no 
simulations of such an event in our sample. 

For the CosmoZo om Ellipticals that have been studied 
in lOser et al.l ()2012l ) we found that the slope of the total 
density correlates with the fraction of stars formed in- 
situ, see Figure [TDJ The more stars have been formed 
within the galaxy itself the steeper the slope of the total 
density, while the accretion of stars by merger events 
flattens the slope. This is in agreement with the fact 
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that for the Binary Ellipticals the steepest slope can be 
found for the 80% gas merger, as discussed above. 

3.3. Evolution of the slopes 

To understand the origin of the total density slopes 
we study the time evolution of the total slopes. In the 
upper panel of Figure [TT] we show the evolution of the 
density and velocity dispersion slopes for all our Binary 
Ellipticals, at the time-steps t = Gyr (initial condition 
time-step), t = 0.58 Gyr, t = 1.75 Gyr, t = 2.32 Gyr and 
t — 2.9 Gyr. Before the merger event occurs we fit the 
slopes to one of the spiral galaxies, in case of an unequal 
mass merger to the more massive spiral. As soon as 
the galaxies are merged we fit the slopes to the remnant 
elliptical. The first passage usually takes place between 
t = 0.5 Gyr and t — 0.6 Gyr, the merger event around 
t = 1.5 Gyr, thus only for the initial condition time t = 
Gyr and the time-step t = 0.58 Gyr we actually fit the 
slopes to the spiral progenitor galaxies. 

Since we set the initial conditions for the Binary Ellip- 
ticals to fit present day spirals, they are already isother- 
mal systems at t — Gyr, i.e., the initial condition slopes 
are around Ttot — —2.0 and /?tot = 0. During the merger 
event, between t = 0.5 Gyr and t = 1.5 Gyr, both the 
total density and velocity dispersion slopes are disturbed 
and not in equilibrium, with flatter slopes for the density 
and slightly positive slopes for the velocity dispersion. 
This shows clearly that the flat density profiles during 
the merger are due to non-equilibrium effects, and that, 
as soon as the central parts of the spirals have merged 
and the central part of the elliptical has formed, both 
slopes return to the isothermal solution and stay (nearly) 
constant. 

Detailed examples for two evolution tracks from t = 
Gyr to t — 3 Gyr are shown in the lower two panels 
of Figure [11] The left panel shows the 1 : 1 spiral merger 
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Fig. 12. — Total velocity slopes versus total density slopes for 
all ellipticals selected from the cosmological re-simulations. Upper 
panel: Slopes at different redshifts for all ellipticals: 2 = 4 (black 
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Lower right panel: Isolated evolution track for the halo M0501-4x, a 
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The black lines in all figures show, like in Figure [3] the analytic 
solution for a spherical system with constant anisotropy. 



with 20% initial gas fraction on a G13 orbit. The right 
panel shows the 3:1 spiral merger with 80% initial gas 
fraction on a G13 orbit, the one elliptical from the binary 
simulations that has a significantly steeper total density 
slope at 3 Gyr than all others. 

The 1:1 spiral merger with 20% initial gas fraction 
shows the typical behavior described above: The merger 
event strongly disturbs the total density slope towards 
natter slopes, while the velocity dispersion slope only 
changes slightly, and after the merger event takes place 
the slopes return to the initial configuration of an isother- 
mal sphere. The merger with an initial gas fraction of 
80% is the only merger event in our sample that be- 
haves differently: It evolves along the black line showing 
the isothermal solutions of the Jeans equation, and the 
slope of the density never reaches equally flat values as 
the merger with 20% initial gas fraction. Even more im- 
portantly, after the merger event takes place the final 
elliptical does not have the same total density and ve- 



locity dispersion slopes as the initial setup. Caused by 
a strong star formation due to gas condensation in the 
center of the newly formed elliptical, as discussed above, 
the final density slope is steeper than the density slope 
of the initial setup. We conclude that a high gas frac- 
tion which causes a significant amount of star formation 
in the central part of the galaxy by condensing gas in 
the galaxies center, is needed for a merger event to cause 
a steeper total density slope once the progenitor system 
has reached an isothermal configuration. 

Figure [T^] shows in the upper panel the evolution of all 
the CosmoZoom Ellipticals, at redshifts z — 4, z = 2, 
z = 0.5, z = 0.2 and z = 0. We always fit the slopes of 
the most massive progenitor of the present day elliptical 
at the given redshift. The lower two panels show the full 
evolution tracks from z — Ato present day for two exam- 
ple ellipticals, one with a multiple merger history (left 
panel) and one with a smooth accretion history, where 
the merger only occur at high redshifts (right panel). 
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We see that the progenitors of the CosmoZoom El- 
lipticals have much steeper density profiles at high red- 
shifts, evolving towards the isothermal spherical system 
case with every merging event, while the velocity dis- 
persion slopes only change slightly during merger events 
and otherwise stay constant. Between z = 4 and z — 2 
the total density slopes do not change much, they are 
all about 7 to t ~ —3, and they do not show an equally 
broad range of values, especially at z = 4, than their 
present day counterparts. After z = 2 the total density 
slopes of all progenitors become flatter and the range of 
slopes broadens. This can be explained by th e two-phase 
forma tion of galaxies that was introduced by Oscr et al. 
(2010): at redshifts z > 2 the formation of the galaxies 
is dominated by the accretion of gas and the in-situ star 
formation followed by a formation phase in which most 
of the mass growth of the galaxies is due to the accretion 
of stars in the form of satellite systems. 

This also explains the correlation between the total 
density slopes and the fraction of stars formed in-situ (see 
Figure ITU]) and the correlation between the total density 
slope and the half-mass radius R\/2 (see Figure [8]). The 
CosmoZoom Ellipticals which have a steep density slope 
have a less dominant mass accretion by mergers in the 
second phase of their formation than those which have 
relatively flat total density slopes, i.e., their in-situ frac- 
tion is higher and they are more compact. CosmoZoom 
Elliptical M0501-4x in the lower right panel of Figure IT21 
is an example for such a merger poor second phase of 
formation. Its total density slope hardly changes from 
z = 4 to present day, and it has no really massive merger 
event in its formation history. 

The opposite case is shown in the lower left panel of 
Figure Q21 CosmoZoom Elliptical M0089-2x has several 
minor merger events with mass ratios between 10:1 and 
3:1 between z = 2 and z = 0.8 and even a major merger 
around z = 0.65. The effect of this major merger event 
can be seen in the lower left panel of Figure [T2] as the 
green circles with the flattest total density slopes. As for 
our Binary Ellipticals, the major merger event disturbs 
the total density slope even beyond 7t t = —2. As for 
the Binary mergers, we also see for the CosmoZoom el- 
lipticals that density slopes lower than 7 to t ~ — 2 are due 
to distortions and non-equilibria during massive merging 
events. 

From both elliptical formation scenarios studied in this 
paper we conclude that merger events, other than in- 
situ star formation, flatten the total density slopes to- 
wards the isothermal solution with a density slope of 
7tot ~ —2. Once the ellipticals have reached this config- 
uration, they stay at this density distribution and, when 
disturbed, evolve back to a density distribution with a 
slope of 7tot ~ — 2, even though the stellar and the dark 
matter component themselves might have changed. We 
therefore conclude that total density distributions with 
a slope of 7tot ~ — 2 act as an attractor solution. Only 
a very gas rich merger event can steepen the slope again 
once it reached the attractor solution, as during a gas 
rich merger the gas can condense in the center of the el- 
liptical and cause significant in-situ star formation, but 
in the second phase of elliptical galaxy formation, which 
is dominated by the accretion and not by the formation 
of stars, a gas rich merger is not a common event. 



4. SUMMARY AND DISCUSSION 

We have investigated a set of 35 spheroidal galaxies 
formed from isolated binary merger events as well as in 
cosmological (zoom) simulation s. The isolated binary 
mer ger spheroids are t aken fr om lJohansson et al.l ((2009a) 
and I Johansson et al.l (|2009b[ ). The spheroids formed in 
cosmological simulations are selected from two different 
simulations: The cosmological zoom-in simul ation sarn- 
ie is a subset of the sample presented in lOser et al] 
2010I) . while the brightest cluster galaxy (BCG) sample 
is selected from a new full hydrodynamical cosmological 
simulation (Dolag et al, in preparation). We analyzed 
the total intrinsic density and velocity dispersion profiles 
of the galaxies to investigate the dependence of the pro- 
files on the formation history, mass, size and dark matter 
fraction. 

We find that all our galaxies are close to isothermal, 
i.e their total (dark matter plus stellar) velocity disper- 
sion distributions are flat, independent of the individ- 
ual total, stellar or dark matter density distributions, 
the mass, the half-mass radius, the dark matter frac- 
tion within one half-mass radius, the environment or the 
formation scenario. The slopes of the total density dis- 
tribution peak at values of Ttot ~ —2.1, with a tendency 
to steeper slopes for less massive, more compact systems 
with lower dark matter fractions within the half-mass 
radius. This is in agreement with obser vational results 
found for the Coma Cluster ellipticals (Thoma s et al.1 
20071 l2009h as well as from strong-lensing (lAueer et al.l 
2O10l:lBarnabe et alJl20ril iSonnenfeld et al.ll2012h ~ 



The cosmological simulations show a similar range of 
values in the density slopes as the observations in agree- 
ment with a similar analysis by iLvskova et al.l (2012) 
for the subset of cosmological zoom in simulations. 
While our isolated binary mergers of present-day galax- 
ies cannot reproduce the observed range of values, bi- 
nary merger between two high-redshift spirals might be 
able to achieve better results since the initial condi- 
tions for those galaxies would look very different. This 
provides evidence to the idea that, at least in a clus- 
ter environment like Coma, the elliptical galaxy pop- 
ulation has not been formed by recent major merg- 
ers but rather at higher redshifts. This is consistent 
with the fact that elliptical galaxies have higher dark 
matter densities than spiral galaxies which indicates a 
higher assembly red shift of about z ~ 3 ( see for example 
iGerhard et al"1l2001t iThbmas et al.ll2009h . Additionally, 
there is strong evidence that massive ellipticals formed 
a signifi cant amount of their stars at high redshifts (for 
example iBrinchmann fc Ellisll2000T ) . supported by obser- 
vatio ns of massive compact ellipticals at high redshift s 
(e.g. Jvan Dokkum et aT1l2009l : lvan de Sande et al.l2011| ). 

There is a small fraction of ellipticals from those 
observations that cannot be reproduced by our simu- 
lated elliptical sample. These objects are character- 
ized by small dark matter fractions, relatively flat to- 
tal density slopes and large dynamical mass-to-light ra- 
tios in the case o f the Coma Cluster ellipticals (see 
iThomas et al.ll20Tll ). A possible explanation for these el- 
lipticals could be that they have a bottom-heavy stellar 
IMF. Non-universal stellar IM Fs have been widely dis - 
cussed recently , for ex a mple by ICappellari et al.l ( 20121) ; 
iFerreras et all (|2012| ): Iron" Dokku m fc Conrovl (|2011l ): 
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IConrov fc van Dokkuml (120121) : Ivan Dokkum fc Conrovl 
<|2012h : iTreu et all (|2010f ). especially in the case of 
massive ellipticals. The predicted dark matter frac- 
tions strongly depend on the assumed IMFs, and as- 
suming for example a Kroupa IMF would result in 
dark matter fractions for the Coma Cluster ellipticals 
that are much higher and thus are in good agreement 
with our simulated results. N eve rtheless, as shown b y 
IConrov fc van Dokkuml (|2012l ) and lWegner et all (|2012fl . 
not even a variable stellar IMF can explain all high dy- 
namical mass-to-light ratios that are observed. 

We also see the effects of adiabatic contraction in 
case of the cosmological simulations: While the den- 
sity slopes for the dark matter component of these el- 
lipticals peak around 7dm ~ —1.67, in good agree- 
ment with the results from ISonncnfcl d et al.1 (2012), 
the dark matter only comparison sample peaks around 
7dm only ~ —1.46. This is con sistent with a num - 
ber of other stud i es, fo r example I Jesseit et all (20021 ) ; 
iBlumenthal et all (119861): iGnedin et all (2004 1201 1[ ) and 
especially I Johansson et al.l ([2012j )~ 

The total density slopes correlate with the fraction of 
stars that are formed in-situ: the steeper the slope, the 
larger the fraction of stars within the galaxy that were 
formed in-situ and the lower the fraction of stars that 
were accreted. This is in agreement with our result that 
the gas fraction in the binary merger scenarios is the 
only component that can significantly alter the slope of 
the total density distribution of the systems. A higher 
star formation rate in the center of the newly formed 
elliptical due to a higher initial gas fraction causes a more 
prominent contribution from the new born stars to the 
total density and thus a steeper total density slope. 

At higher redshifts, where gas and in-situ star forma- 
tion dominate the galaxies, the ellipticals from cosmolog- 
ical simulations have a total density slope of jtot ~ — 3, 
evolving through merger events towards a slope of jtot ~ 
—2, supp orting the idea o f the two-phase formation of 
galaxies (|Oser et al.l l2010f) . If the in-situ fraction of a 
galaxy at z—0 is high, then the galaxy has accreted less 
stars from its environment in the second phase of its for- 
mation than a galaxy with a low in-situ fraction at z=0. 
Without enough gas-poor accretion in the second phase 
of formation, the total density slope could not changed a 
lot towards 7 to t = —2, and thus the slope stays close to 
7t ot ~ —3. This is in agr eement with results presented 
in Uohansson et al.l (20121 ) but in disagreement wit h ob- 
servations from stro ng lensing bv iRuff et al.l ()2011l ) and 
iBolton et al.l (2012ft who find a slight trend indicating 
that the slopes of the total density profiles of early-type 
galaxies at higher redshift are slightly flatter. Never- 
theless, this observed trend is very weak, and it could 
be due to a bias in the weak lensing sample to wards 
merging systems, as reported bv lTorri et~a l. (2004), who 
found that merging systems tend to boost strong lens- 
ing. An enhancement in the lensing efficiency reported 



by IZitrin et al.l (2012ft for an observed merging cluster of 
galaxies also supports this idea. This would be in agree- 
ment with our result that merging systems show much 
flatter density slopes. 

Our simple model predicts that the steepness of the 
slope of present day galaxies is a measurement for the 
importance of the mergers the elliptical galaxy went 
through in its second formation phase. Ellipticals with 
steep slopes close to 7tot ~ — 3 had (nearly) no merger 
event in this second phase, while ellipticals with slopes 
around 7tot ~ — 2 had a strong, collisionless merger dom- 
inated second formation phase. Since our ellipticals with 
slopes around 7tot ~ — 2 have generally higher dark mat- 
ter fi^£tiojis, tlris_isconsistent with the results presented 
bv iHilz et al.l (2012b. a) who show that the accretion of 
several small satellite systems strongly increases the dark 
matter fraction within the half-mass radius. From both, 
binary merger ellipticals and ellipticals from cosmologi- 
cal framework, we see that, once an elliptical has reached 
a total density slope of jtot ~ — 2, further merger events 
do not change the slope anymore. 

We conclude that the density distributions with a slope 
of 7tot ~ — 2 acts as attractors, independent of the indi- 
vidual stellar mass distributions of ellipticals. We suggest 
that all elliptical galaxies will in time end up in such a 
configuration. However, the mechanism that leads to this 
attractor is still unclear. A possible explanation could be 
that this attractor state is a result of violent re laxation 
(Lvnden-Belllll967ft . although |Snz_eLalJ (2012bft showed 
that violent relaxation is less efficient in minor merger 
events then in major merger events. 

In general, the relaxation times for the elliptical galax- 
ies from cosmological frameworks are much faster then 
for isolated binary mergers due to the presence of sub- 
structures and potential fluctuations. We conclude that, 
to understand the entire range of elliptical galaxies ob- 
served at present day and their complex evolution sce- 
narios, a full cosmological treatment is needed. 
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